Entrance Channel Dynamics of Hot and Cold Fusion Reactions 
Leading to Superheavy Elements 
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We investigate the entrance channel dynamics for the reactions ^"Zn + ^"^^Pb and **Ca + ^-^^U 
using the fully microscopic time-dependent Hartree-Fock (TDHF) theory coupled with a density 
constraint. We calculate excitation energies and capture cross-sections relevant for the study of 
superheavy formations. We discuss the deformation dependence of the ion-ion potential for the 
**Ca -I- 2'^*U system and perform an alignment angle averaging for the calculation of the capture 
cross-section. The results show that this parameter-free approach can generate results in good 
agreement with experiment and other theories. 
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I. INTRODUCTION 



sections than those predicted by zero temperature cal- 
culations [12)]. 
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One of the most fascinating research areas involv- 
ing low-energy nuclear reactions is the search for super- 
heavy elements. Experimentally, two approaches have 
been used for the synthesis of these elements, one uti- 
lizing closed shell nuclei with lead-based targets (cold- 
fusion) [i|,i2, the other utilizing deformed actinide tar- 
gets with "'^Ca projectiles (hot-fusion) While both 
methods have been successful in synthesizing new ele- 
ments the evaporation residue cross-sections of the hot- 
fusion reactions were found to be as much as three times 
larger than those of the cold fusion ones. To pinpoint 
the root of this difference it is important to understand 
the details of the entrance channel dynamics of these 
systems since the properties of the dinuclear system at 
the capture point will strongly influence the outcome of 
the reaction. For light and medium mass systems the 
capture cross-section may be considered to be the same 
as that for complete fusion, whereas for heavy systems 
leading to superheavy formations the evaporation residue 
cross-section is dramatically reduced due to the quasi- 
fission and fusion-fission processes, thus making the cap- 
ture cross-section to be essentially the sum of these two 
cross-sections. What is also difficult to ascertain is the 
configuration of the composite system namely, whether 
the system has a single-center compound-like configura- 
tion or a dinuclear configuration acconipanied by parti- 
cle exchange. Most dynamical models |6l-[l]| argue that 
for heavy systems a dinuclear complex is formed initially 
and the barrier structure and the excitation energy of 
this precompound system will determine its survival to 
breaking up via quasi-fission. Furthermore, if the nucleus 
survives this initial state and evolves to a compound sys- 
tem it can still fission due to its excitation. Recent mi- 
croscopic calculations have shown that the temperature 
(excitation energy) of the compound systems strongly al- 
ters the barrier structures of potential energy surfaces, 
with implications of very different fusion-fission cross- 



It is generally acknowledged that the TDHF theory 
provides a useful foundation for a fully microscopic many- 
body theory of low-energy heavy-ion reactions (l3| . Re- 
cent three-dimensional TDHF calculations with no sym- 
metry assumptions and using modern Skyrme forces have 
been shown to accurately reproduce phenomena deter- 
mined by the early stages of the heavy-ion dynamics [ij- 
ITgI I . Recently, we have developed the density-constrained 
TDHF (DC-TDHF) method which is based on the 
generalization of the density constraint method devel- 
oped earlier We have shown that using the DC- 
TDHF method ion-ion potential barriers can be accu- 
rately produced in most cases (l9l - l2l| . as these calcula- 
tions also depend on early stages of the ion-ion dynamics. 
Furthermore, one-body energy dissipation extracted from 
TDHF for low-energy fusion reactions was found to be in 
agreement with the friction coefficients based on the lin- 
ear response theory as well as those in models where the 
dissipation was specifically adjusted to describe experi- 
ments [23] . All of these new results suggest that TDHF 
dynamics may provide a good description of the early 
stages of heavy-ion collisions. 



Recently, we have also introduced an approach to ex- 
tract excitation energies directly from full microscopic 
TDHF calculations [23|. In this manuscript we perform 
TDHF calculations accompanied by density constraint 
calculations for ^°Zn -I- ^"*Pb and ''^Ca -I- ^^*^U systems, 
which represent typical examples of cold and hot fu- 
sion reactions leading to superheavy formation, respec- 
tively. In addition to calculating the excitation energy at 
the capture point we also investigate the capture cross- 
section and try to elucidate the differences between these 
two reactions. 
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II. THEORETICAL OUTLINE 

In this Section we will discuss issues pertaining to the 
application of TDHF to study collisions involving heavy 
reaction partners. 



A. TDHF Dynamics 

It is generally accepted that TDHF theory is a candi- 
date for a microscopic theory that may provide a unified 
approach to the description of diverse physical phenom- 
ena such as fusion, deep-inelastic collisions, dinuclear and 
compound nucleus formation, and possibly fission. Since 
TDHF is based on the independent-particle approxima- 
tion it can be interpreted as the semi-classical limit of a 
fully quantal theory thus allowing a connection to macro- 
scopic coordinates and providing insight about the colli- 
sion process. In this sense the TDHF dynamics can only 
compute the semiclassical trajectories of the collective 
moments of the composite system as a function of time. 
The presence of residual interaction, absent in TDHF, 
may produce fluctuations and correlations which effect 
the mean values of these trajectories. 

For TDHF collisions of light and medium mass sys- 
tems as well as highly mass-asymmetric systems fusion 
generally occurs immediately above the Coulomb bar- 
rier, while in heavier systems there is an energy range 
above the barrier where fusion does not occur. This phe- 
nomenon is the microscopic analogue of the macroscopic 
extra-push threshold (23 | and has been recently studied 
for TDHF collisions involving heavy and nearly symmet- 
ric reaction partners [25| . In the lower part of this en- 
ergy range deep-inelastic collisions are dominant while at 
slightly higher energies the system develops a long-lived 
and pronounced neck reminiscent of a dinuclear configu- 
ration. The outcome of these long-lived configurations is 
uncertain due to the absence of quantal shell effects and 
fluctuations. For these systems TDHF results in a com- 
pact configuration only for energies considerably above 
the static potential energy surface. However, despite the 
high energy, single-particle friction can quickly absorb 
this energy and lead to a configuration that may be con- 
sidered a thermal doorway state. As long as the average 
single-particle excitation energy per nucleon in this door- 
way state is less than the shell energy (about 4 — 8 MeV) 
the details of the ground state potential energy surface 
are still felt and shell correction energies infiuence the 
TDHF dynamics. It is precisely for this reason that the 
DC-TDHF approach allows us to reproduce ion-ion in- 
teraction barriers for heavy-ion collisions. 



B. DC-TDHF and Excitation Energy 



used to perform a static Hartree-Fock minimization while 
holding the neutron and proton densities constrained to 
be the corresponding instantaneous TDHF densities. In 
essence, this provides us with the TDHF dynamical path 
in relation to the multi-dimensional static energy sur- 
face of the combined nuclear system. The advantages of 
this method in comparison to other mean-field based mi- 
croscopic methods such as the constrained Hartree-Fock 
(CHF) method are obvious. First, there is no need to 
introduce artificial constraining operators which assume 
that the collective motion is confined to the constrained 
phase space: second, the static adiabatic approximation 
is replaced by the dynamical analogue where the most 
energetically favorable state is obtained by including sud- 
den rearrangements and the dynamical system does not 
have to move along the valley of the potential energy 
surface. In short we have a self-organizing system which 
selects its evolutionary path by itself following the micro- 
scopic dynamics. All of the dynamical features included 
in TDHF are naturally included in the DC-TDHF cal- 
culations. These effects include neck formation, mass 
exchange, internal excitations, deformation effects to all 
order, as well as the effect of nuclear alignment for de- 
formed systems. In the DC-TDHF method the ion-ion 
interaction potential is given by 



ViR) = Euc{R)-EA,-EA, 



(1) 



In the DC-TDHF approach [l7l the TDHF time- 
evolution takes place with no restrictions. At certain 
times during the evolution the instantaneous density is 



where -Edc is the density-constrained energy at the in- 
stantaneous separation R(t), while _Bai and are the 
binding energies of the two nuclei obtained with the same 
effective interaction. In writing Eq. ([T} we have intro- 
duced the concept of an adiabatic reference state for a 
given TDHF state. The difference between these two en- 
ergies represents the internal energy. The adiabatic refer- 
ence state is the one obtained via the density constraint 
calculation, which is the Slater determinant with lowest 
energy for the given density with vanishing current and 
approximates the collective potential energy 18]. We 
would like to emphasize again that this procedure does 
not affect the TDHF time-evolution and contains no free 
parameters or normalization. 

Ion-ion interaction potentials calculated using DC- 
TDHF correspond to the configuration attained during a 
particular TDHF collision. For light and medium mass 
systems as well as heavier systems for which fusion is the 
dominant reaction product, DC-TDHF gives the fusion 
barrier with an appreciable but relatively small energy 
dependence. On the other hand, for reactions leading 
to superheavy systems fusion is not the dominant chan- 
nel at barrier top energies. Instead the system sticks in 
some dinuclear configuration with possible break-up af- 
ter exchanging a few nucleons. The long-time evolution 
to break-up is beyond the scope of TDHF due to the 
absence of quantal shell effects and fluctuations. As we 
increase the energy above the barrier this phenomenon 
gradually changes to the formation of a truly compos- 
ite object. This is somewhat similar to the extra-push 
phenomenon discussed in phenomenological models. For 
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this reason the energy dependence of the DC-TDHF in- 
teraction barriers for these systems is not just due to the 
dynamical effects for the same final configuration but ac- 
tually represent different final configurations. 

The calculation of the excitation energy is achieved by 
dividing the TDHF motion into a collective and intrinsic 
part [23|. The major assumption in achieving this goal 
is to assume that the collective part is primarily deter- 
mined by the density p(r, t) and the current j(r, t). Con- 
sequently, the excitation energy can be formally written 
as 

E*it) - Etdhf - Ecoii ip{t),m) , (2) 

where -Etdhf is the total energy of the dynamical system, 
which is a conserved quantity, and E'coii represents the 
collective energy of the system. In the next step we break 
up the collective energy into two parts 

Ecou (t) = Ek.,n ipit),m) + Edc (pit)) , (3) 
where i?kin represents the kinetic part and is given by 

i?H„(pWJW) = y J d\ fit)/ Pit), (4) 

which is asymptotically equivalent to the kinetic energy 
of the relative motion, ^fiR^, where p, is the reduced 
mass and is the ion-ion separation distance. The 
dynamics of the ion-ion separation is provided by an 
unrestricted TDHF run thus allowing us to deduce the 
excitation energy as a function of the distance parameter, 
E*iR). 

III. RESULTS 

Calculations were done in 3-D geometry and using the 
full Skyrme force (SLy4) [2^ without the centcr-of-mass 
correction as described in Ref. [27j. We have performed 
density constraint calculations every 20 time steps. For 
the calculation of the ion-ion separation distance R we 
use the hybrid method, which relates the coordinate to 
the quadrupole moment for small R values, as described 
in Ref. (2i|. The accuracy of the density constraint cal- 
culations is commensurate with the accuracy of the static 
calculations. 



A. ^'^Ca+^^'^u System 

As an example of superheavy formation from a hot- 
fusion reaction we have studied the '^^Ca-l-^^^U system. 
Hartree-Fock (HF) calculations produce a spherical '^^Ca 
nucleus, whereas ^^Sy j^g^g large axial deformation. The 
large deformation of ^ssy jg expected to strongly influ- 
ence the interaction barriers for this system. This is 
shown in Fig. [1] which shows the interaction barriers, 
ViR), calculated using the DC-TDHF method as a func- 
tion of cm. energy and for three different orientations 



of the 238 u nucleus. The alignment angle /3 is the angle 
between the symmetry axis of the ^38^ nucleus and the 
collision axis. Also shown in Fig. [1] is the point Coulomb 
potential corresponding to this collision. The deviations 
from the point Coulomb potential at large R values are 
due to the deformation of the ^38^ nucleus. We first 
notice that the barriers corresponding to the polar ori- 
entation i(3 = 0°) of the 238^ nucleus are much lower 
and peak at larger ion-ion separation distance R. On the 
other hand, the barriers corresponding to the equatorial 
orientation of ^38^ (yj = 90°) are much higher and peak 
at smaller R values. For the intermediate values of (3 the 
barriers rise rapidly as we increase the orientation angle 
from /3 — 0°, as can be seen for (3 = 45°. The rise in the 
barrier height as a function of increasing /3 values is not 
linear but seems to rise more rapidly for smaller /3 val- 
ues. We also see that for lower energies central collisions 
with polar orientation of ^38^ ^^.g ^]-^g only orientations 
which result in the sticking of the two nuclei, while the 
equatorial orientations of ^38^ result in a deep-inelastic 
collision. Also, shown in Fig. [T] are the experimental en- 
ergies d, Q for this reaction. We observe that all of the 
experimental energies are above the barriers obtained for 
the polar alignment of the ^38^ nucleus. 




R(fm) 

FIG. 1: (color online) Potential barriers, V{R), for the 
■'^'Ca-t-^^^U system obtained from DC-TDHF calculations us- 
ing Eq. IT} as a function of Ec.m. energy and for selected 
orientation angles (3 of the ^^®U nucleus. Also, shown are the 
experimental cm. energies. 

Furthermore, the potentials shown in Fig. [T] display a 
very strong energy dependence. Detailed analysis of the 
TDHF time-evolution and density profiles show that at 
lower cm. energies and for central collisions the polar 
configuration of ^38^ leads to a dinuclear system where 
both nuclei maintain their cores and exchange nucleons. 
Non-central collisions at these energies result in deep- 
inelastic fragments. As we mentioned above, for these low 
energy collisions the equatorial collisions result in deep- 
inelastic reaction products even for central collisions. At 
higher energies the system forms a true composite with 
overlapping cores. These are the potentials which should 
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be used for calculating capture cross-sections. 

For the formation of a superheavy system the ex- 
citation energy carries great significance since a high- 
excitation energy at the capture point would result in 
quasi-fission events while high excitation of the com- 
pound nucleus will lead to fusion-fission. Many factors 
play a role in building up the excitation energy. As we 
have discussed in the Introduction, recent TDHF calcu- 
lations suggest that during the early phase of the ion-ion 
collisions TDHF theory may provide a good approxima- 
tion for the transfer of the initial kinetic energy into inter- 
nal and collective excitations via the dynamical evolution 
of shell effects. Naturally, the excitation also depends on 
other nuclear properties such as deformation and align- 
ment, mass asymmetry in the entrance channel, and im- 
pact parameter. In Fig. [2] we show the excitation energy 
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FIG. 2: (color online) Excitation energy, E*{R), in a central 
collision for various -Ec.m. energies and for alignment angles 
angles /3 = 0° and /3 = 90° of the ^^*U nucleus. 



E* (R) as a function of cm. energy and for two alignment 
angles (/3 = 0° and ;3 = 90° ) of the ^^^\J nucleus. The 
excitation energy curves start at zero excitation when 
the two nuclei are far apart, which also provides a test 
for the numerical accuracy of the calculation. We note 
that the system is excited much earlier during the colli- 
sion process for the polar alignment of the ^ssu nucleus 
and has a higher excitation than the corresponding col- 
lision for the equatorial orientation. Only two curves are 
shown for the equatorial collision since at lower energies 
we have deep- inelastic collisions for this alignment. We 
note that the highest point reached for these excitation 
curves is chosen to be the point where the nuclei almost 
come to a stop inside the barrier, which corresponds to 
a nearly zero collective kinetic energy. For energies for 
which the collision outcome is capture this would be the 
excitation energy at the capture point. 

In Fig. [3] we show the excitation energy as a func- 
tion of the center-of-mass energy for both ^^Ca-I-^^^U 
and ™Zn-|-^°*Pb systems calculated at the potential min- 
imum inside the barrier. The excitation energy curves in- 
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FIG. 3: (color online) Excitation energy in a central collision 
at the potential minimum inside the barrier, E*, as a function 
of iJc.m. for the '^^Ca+^^'^U and the ™Zn+2''**Pb systems. For 
the case of ^*Ca-|-^'^*U system ^•^^U alignment angles of /? = 
0°, 45°, and (3 = 90° are shown. 



crease with a typical linear slope as a function of Ec.m.- 
For the case of ^^Ca-I-^'^^U system we show the excita- 
tion energies for ^ssy alignment angles of /3 = 0°,45°, 
and /3 = 90°. The upward arrows in the bottom of the 
figure denote the experimental energies for these reac- 
tions. For larger alignment angles of ^■^^U, the curves 
start at a higher center-of-mass energy since at lower en- 
ergies we only have deep-inelastic collisions. The lowest 
two points for the ™Zn-|-^°^Pb curve (shown as stars) are 
simply the extrapolated values down to the experimental 
energies. 

While at first glance the excitation energies for the 
^^Ca-f^'^^U system look much higher than those for the 
™Zn-|-^°^Pb system, this may be somewhat misleading. 
The reason is that the excitations for the ^^Ca-I-^'^^U 
system must be angle averaged for different alignments 
of the 238u nucleus and also weighted by the alignment 
probability discussed below. So, for example the (3 = 0° 
excitation at a particular center-of-mass energy will be 
multiplied by the sin(/?) factor in the integration weight 
and will make a zero contribution. The value of sin(/3) in- 
creases with l3, but at the same time the excitation energy 
decreases. Since TDHF calculations at the experimental 
energies only yield dinuclear configurations for the small 
values of /3 but rapidly move into the deep-inelastic do- 
main for larger values we could not compute this average 
in practice. Naturally, a fully quantal system will have 
a certain probability for resulting in a dinuclear config- 
uration for all values of sm{/3) at these energies and the 
excitation energy will be smaller than the one one shown 
for the (3 = 0° case. 

As we have discussed earlier, for systems leading to su- 
perheavy formation the evaporation residue cross-section 
is customarily represented in terms of the various phases 
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of the reaction process as 



O'ER 



-'capture 



PcN • R 



survival 



(5) 



where cter denotes the evaporation residue cross-section 
for the superheavy system, cTcapturo is the capture cross- 
section for the two-ions, Pcn is the probabihty of form- 
ing a compound nucleus, and Pgurvivai is the probabil- 
ity that this compound system survives various breakup 
and fission events. The calculations presented here can 
only address the capture cross-section for these systems 
since the subsequent reaction possibilities are beyond the 
scope of the TDHF theory. For most light systems for 
which fusion is the dominant reaction product (Tcapturo 
and O'ER are essentially the same and equal to the fusion 
cross-section, CTfusion- Instead, for reactions involving su- 
perheavy formations we have 



''capture — O-QF + CTpF + O'ER , 



(6) 



where ctqf and uff denote the quasi-fission and fusion- 
fission cross-sections, respectively. For these reactions 
the evaporation residue cross-section, (Ter, is very small 
and therefor the capture cross-section is to a large ex- 
tent equal to sum of the two fission cross-sections. Fur- 
thermore, the distinction between deep-inelastic reac- 
tions and quasi-fission is somewhat difficult and usu- 
ally achieved by setting windows for fragment masses of 
Af = Acn/2 ± 20 and on their kinetic energy. 



where dP/(d/3 sin(/3)) represents the alignment proba- 
bility and (j{Ec.m.iP) is the capture cross-section asso- 
ciated with a particular alignment. In practice, we cal- 
culate the alignment probability due to Coulomb excita- 
tion as described in Ref. 31]. Due to the relatively small 
charge of the *®Ca nucleus this probability is in the range 
0.46 — 0.52 as shown in Fig. Uand does not vary appre- 
ciably with energy. One important fact to notice in the 
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FIG. 5: 

''*Ca-|-^^*U system as a function of Ec. 



(color online) Capture cross-sections for the 

energy (black cir- 
cles). Also, shown are the experimental cross-sections (red 
squares) [J. 
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FIG. 4: Dynamic alignment due to Coulomb excitation of 
^^*U. Shown is the orientation probability as a function of the 
Euler angle /3 in a central collision at internuclear distances 
R = 1500 fm (dashed horizontal curve) and at _R = 20 fm at 
center-of-mass energy of i?c.m. = 220 MeV (solid curve). 



cross-section formula given in Eq. ([T]) is that the cross- 
section is multiplied by the sin(/3) factor, which renders 
the contribution originating from the lowest barriers at 
small values of (3 to be very small. However, unlike the 
calculation of excitation energy, which requires a TDHF 
collision at exactly the same experimental energy, the po- 
tential barriers are obtained by performing TDHF calcu- 
lations at higher energies and using the density constraint 
to calculate the ion-ion potentials. For a consistent calcu- 
lation of fusion cross sections above and below the barrier 
energies we have adopted the commonly used Incoming 
Wave Boundary Condition (IWBC) method [H,!!!]. In 
practice, we have varied the alignment angle /3 in 10° 
steps between 0° and 90°. In Fig. [S] we show the capture 
cross-sections for the ^^Ca-|-^'^*U system as a function of 
^-c.m. energy (black circles). Also, shown are the experi- 
mental cross-sections (red squares) [^] . While the agree- 
ment with experiment is not perfect it is still remarkably 
good since these calculations contain no free parameters. 



For the calculation of the capture cross-section we need 
to av erage o ver all possible alignments of the ^ssy 
cleus |28l - {33 | as given by 



capture 



d/3 sin(/3) ;t/^7^ o{E, 
d^ sm(/3) 



,fi) , 
(7) 



B. ™Zn+2'*Pb System 

We have also investigated the cold-fusion reaction of 
™Zn-|-^°^Pb, which lead to the same superheavy ele- 
ment Z = 112 (but a different isotope). The Hartree- 
Fock calculations for ^Zn and ^"^Pb produce a spheri- 
cal nucleus in both cases. In Fig. [B] we plot the potcn- 
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tial barriers for the ^"Zn+^^^Pb system obtained from 
DC-TDHF calculations using Eq. ([T]) as a function of 
Ec.m. energy. Also shown is the point Coulomb potential 
for two spherical nuclei. Unlike light nuclei, the devia- 
tions from the Coulomb trajectory start relatively early 
(15 — 16 fm) due to the large charge of the two ions. A 
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FIG. 6; (color online) Potential barriers, V{R), for the 
™Zn+^"*Pb system obtained from DC-TDHF calculations us- 
ing Eq. (O, for various -Ec.m. energies as indicated. 

similar behavior to the ^*Ca-t-^^^U system is observed for 
the ^"Zn-t-^^^Pb system, namely at lower energies nuclear 
densities show a dinuclear character for central collisions 
and result in deep-inelastic collision at non-central im- 
pact parameters. The lowest barrier in Fig. [S] shows the 
case at Ec.m. = 265 MeV where even for a central collision 
the result is deep-inelastic as can be seen from the dotted 
blue curve. Only at the highest two energies the densities 
corresponding to central collisions show a more composite 
character, which may be identified as capture, and result 
in a similar density distribution for non-central collisions 
as well. We also observe the flattening of the ion-ion 
potentials as the nuclear overlap increases at smaller R 
values, which seems to be a general behavior for collisions 
of two heavy nuclei. In Fig. [7] we show the excitation en- 
ergy, E*{R), as a function of the E'c.m. energy for the 
^°Zn-t-^°^Pb system. Again, the highest value attained 
is at the point where the two nuclei come to a stop in- 
side the ion-ion potential. These energies were shown in 
Fig. [3] as a function of Ec.m. energy. 

We have also calculated the capture cross-section for 
the ™Zn-|-^°^Pb system. Since we are dealing with two 
spherical nuclei no angle averaging is necessary for this 
system. However, we could not find experimental mea- 
surements for the capture cross-section. For comparison 
we have used a model calculation where the cross-sections 
for the synthesis of superheavy elements were analyzed 
using the concept of a dinuclear system [2^ . The authors 
calculated capture or quasi-fission, fusion, and In evap- 
oration residue cross-sections which reproduced a single 
In evaporation residue experimental cross-section. For 
comparison we have plotted our capture cross-sections 
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FIG. 7: (color online) Excitation energy, E*{R) for the 
^"Zn-l-^'^^Pb system and for various -Ec.m. energies as indi- 
cated. 



together with their findings for the two experimental lab- 
oratory energies ^ of -Ejab = 343.8 MeV and 346.1 MeV 
in Fig. [51 Our calculations at the two experimental en- 
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FIG. 8: (color online) Cross-sections (in ^h) from Ref. [2y] 
as a function of -Biab energy for the ^"Zn-t-^^^Pb system. The 
dotted curve shows the quasi-fission, the dashed curve the 
fusion, and the dot-dashed curve the In evaporation residue 
cross-section. Our calculations at the two experimental en- 
ergies corresponding to the two potential energies leading to 
capture in Fig. [6] are shown as filled circles. 

ergies corresponding to the two potential energy curves 
leading to capture in Fig. [B] are shown as filled circles. 
Differences in laboratory energies between our results 
and the results of the model calculation f29j stem from 
the fact that model calculations make use of the com- 
pound nucleus excitation energy and relate this to the 
laboratory energy. We have directly used the experimen- 
tal energies from Ref. since we do not compute the 
compound nucleus excitation energy. On the other hand, 
the dependence of the quasi-fission cross-section on lab- 
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oratory energy is very flat as shown in Fig. [8] As we 
see the calculated cross-sections at the higher laboratory 
energy of 346.1 MeV is in a reasonable agreement with 
the model calculation for both potential energy curves, 
whereas for the lower laboratory energy of 343.8 MeV 
the lower potential energy curve substantially underes- 
timates the model calculation. We may conclude from 
this observation that the higher potential energy curve 
provides a better representation of the inner part of the 
ion-ion barrier. However, since we cannot rule out the 
lower potential energy curve as a capture event based on 
our TDHF results we have shown both results. 



IV. SUMMARY 

In this work we have investigated two systems which 
are known to produce the element Z = 112 in experi- 
ments. We find that the collisions of such heavy systems 
have very different characteristics than the TDHF calcu- 
lations of light and medium-mass systems. We have used 
the density constraint along with TDHF to obtain ion- 
ion interaction potentials and excitation energies. The 
dependence of the ion-ion potential on the deformation 
of the ^^®U nucleus was studied. The calculated capture 



cross-sections are found to be in reasonable agreement 
with data and other model calculations. 

The fully microscopic TDHF theory has shown it- 
self to be rich in nuclear phenomena and continues to 
stimulate our understanding of nuclear dynamics. The 
time-dependent mean-field studies seem to show that the 
dynamic evolution builds up correlations that are not 
present in the static theory. While modern Skyrme forces 
provide a much better description of static nuclear prop- 
erties in comparison to the earlier parametrizations there 
is a need to obtain even better parametrizations that in- 
corporate deformation and scattering data into the fit 
process. 
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